arXiv:1501.04395vl [cs.IT] 19 Jan 2015 


IEEE TRANSACTIONS ON VEHICULAR TECHNOLOGY SUBMISSION 


1 


Spherical Harmonic Expansion of Fisher-Bingham 
Distribution and 3D Spatial Fading Correlation for 

Multiple-Antenna Systems 

Yibeltal F. Alem, Student Member, IEEE, Zubair Khalid, Member, IEEE and Rodney A. Kennedy, Fellow, IEEE 


Abstract —This paper considers the 3D spatial fading correla¬ 
tion (SFC) resulting from an angle-of-arrival (AoA) distribution 
that can be modelled by a mixture of Fisher-Bingham distribu¬ 
tions on the sphere. By deriving a closed-form expression for the 
spherical harmonic transform for the component Fisher-Bingham 
distributions, with arbitrary parameter values, we obtain a 
closed-form expression of the 3D-SFC for the mixture case. The 
3D-SFC expression is general and can be used in arbitrary multi¬ 
antenna array geometries and is demonstrated for the cases of 
a 2D uniform circular array in the horizontal plane and a 3D 
regular dodecahedral array. In computational aspects, we use 
recursions to compute the spherical harmonic coefficients and 
give pragmatic guidelines on the truncation size in the series 
representations to yield machine precision accuracy results. The 
results are further corroborated through numerical experiments 
to demonstrate that the closed-form expressions yield the same 
results as significantly more computationally expensive numerical 
integration methods. 

Index Terms —Fisher-Bingham distribution, spherical har¬ 
monic expansion, spatial correlation, MIMO, angle of ar¬ 
rival (AoA). 

I. Introduction 

The Fisher-Bingham distribution, also known as the Kent 
distribution, belongs to the family of spherical distributions 
in directional statistics jT|. It has been used in applications 
in a wide range of disciplines for modelling and analysing 
directional data. These applications include sound source lo¬ 
calization |2j, |3j, joint set identification |4j, modelling protein 
structures {5), 3D beamforming |6|, classification of remote 
sensing data 0, modelling the distribution of AoA in wireless 
communication 0, to name a few. 

In this work, we focus our attention to the use of Fisher- 
Bingham distribution for modelling the distribution of angle 
of arrival (AoA) (also referred as the distribution of scatterers) 
in wireless communication and computation of spatial fading 
correlation (SFC) experienced between elements of multiple- 
antenna array systems. Modelling the distribution of scatterers 
and characterising the spatial correlation of fading channels 
is a key factor in evaluating the performance of wireless 
communication systems with multiple antenna elements fT0|- 
G3- It has been an active area of research for the past two 
decades or so and a number of spatial correlation models and 
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closed-form expressions for evaluating the SFC function have 
been developed in the existing literature (e.g., (TO), (1§-J25)). 

The elliptic (directional) nature of the Fisher-Bingham 
distribution offers flexibility in modelling the distribution of 
normalized power or AoA of the multipath components for 
most practical scenarios. Exploiting this fact, a 3D spatial 
correlation model has been developed in (9), where the dis¬ 
tribution of AoA of the multipath components is modelled by 
a positive linear sum of Fisher-Bingham distributions, each 
with different parameters. Such modelling has been shown 
to be useful in a sense that it fits well with the multi-input 
multi-output (MIMO) field data and allows the evaluation 
of the SFC as a function of the angular spread, ovalness 
parameter, azimuth, elevation, and prior contribution of each 
cluster. Although the SFC function presented in |9| is general 
in a sense that it is valid for any arbitrary antenna array 
geometry, it has not been expressed in closed-form and has 
been only evaluated using numerical integration techniques, 
which can be computationally intensive to produce sufficiently 
accurate results. If the spherical harmonic expansion of the 
Fisher-Bingham distribution is given in a closed-form, the 
SFC function can be computed analytically using the spherical 
harmonic expansion of the distribution of AoA of the multipath 
components To the best of our knowledge, the spherical 
harmonic expansion of Fisher-Bingham distribution has not 
been derived in the existing literature. 

In the current work, we present the spherical harmonic 
expansion of Fisher-Bingham distribution and a closed-form 
expression that enables the analytic computation of the spher¬ 
ical harmonic coefficients. We also address the computational 
considerations required to be taken into account in the evalua¬ 
tion of the proposed closed-form. Using the proposed spherical 
harmonic expansion of the Fisher-Bingham distribution, we 
also formulate the SFC function experienced between two 
arbitrary points in 3D-space for the case when the distribution 
of AoA of the multipath components is modelled by a mixture 
(positive linear sum) of Fisher-Bingham distributions. The 
SFC presented here is general in a sense that it is expressed as 
a function of arbitrary points in 3D-space and therefore can be 
used to compute spatial correlation for any 2D and 3D antenna 
array geometries. Through numerical analysis, we also validate 
the correctness of the proposed spherical harmonic expansion 
of Fisher-Bingham distribution and the SFC function. In this 
paper, our main objective is to employ the proposed spherical 
harmonic expansion of the Fisher-Bingham distribution for 
computing the spatial correlation. However, we expect that 
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the proposed spherical harmonic expansion can be useful in 
various applications where the Fisher-Bingham distribution is 
used to model and analyse directional data (e.g., @-0). 

The remainder of the paper is structured as follows. We 
review the mathematical background related to signals defined 
on the 2-sphere and spherical harmonics in Section [II] In 
Section [HI] we define the Fisher-Bingham distribution and 
present its application in modelling 3D spatial correlation. 
We present the spherical harmonic expansion of the Fisher- 
Bingham distribution and analytical formula for computing 


NT = 


4-7T (l + m)V 


' (£ + m)\ 
(£ — m)\ 


l m, o 


W, 


where d ^ m , (•) denotes the Wigner-d function of degree £ and 
orders m and m! [ 261. 

By completeness of spherical harmonics, any finite energy 
function f{x) on the 2-sphere can be expanded as 


/(£) = £ £(/)r*7*(£). 


(4) 


£—0 m——£ 


the spherical harmonic coefficients in Section IV where we 
also address computational issues. We derive a closed-form 
expression for the 3D SFC between two arbitrary points in 30- 
space when the AoA of an incident signal follows the Fisher- 
Bingham probability density function (pdf) in Section [V] 
In Section |VI| we carry out a numerical validation of the 
proposed results and provide examples of the SFC for uniform 
circular array (2D) and regular dodecahedron array (3D) 
antenna elements. Finally, the concluding remarks are made 
in Section [VII] 

II. Mathematical Preliminaries 

A. Signals on the 2-Sphere 

We consider complex valued square integrable functions 
defined on the 2-sphere, 8 2 . The set of such functions forms 
a Hilbert space, denoted by L 2 (§ 2 ), that is is equipped with 
the inner product defined for two functions / and g defined 

on S 2 as (261 

{f,9)= [ f{x)g{x)ds{x ), (1) 

J s 2 

which induces a norm ||/|| = (/, f) 1 ^. Here, (•) 

denotes the complex conjugate operation and x = 
[sin0cos0,sin0sin0,cos#] 7 £ § 2 C R 3 represents a point 
on the 2-sphere, where [•] 7 represents the vector transpose, 
9 € [0,7r] and £ [0, 27t) denote the co-latitude and longitude 
respectively and ds(i r) = sin 0 dO dip is the surface measure on 
the 2-sphere. The functions with finite energy (induced norm) 
are referred as signals on the sphere. 

B. Spherical Harmonics 

Spherical harmonics serve as orthonormal basis functions 
for the representation of functions on the sphere and are 
defined for integer degree £ > 0 and integer order \rn\ < £ as 

Y<r{x) = Y™(9, fi) = N™P ™(cos 9)e im ^, 

with 

A he + i (£ — m)! 


where (/)” is the spherical harmonic coefficient given by 

{/)?=[ f{x)Y e m (x)ds(x). (5) 

J s 2 

The signal / is said to be band-limited in the spectral domain 
at degree L if (/)™ = 0 for t > L. For a real-valued function 
/(a;), we also note the relation 


(/)7 m = (-1 r(f)T, 


( 6 ) 


which stems from the conjugate symmetry property of spher¬ 
ical harmonics 


C. Rotation on the Sphere 

The rotation group SO(3) is characterized by Euler angles 
(<p,d, to) £ SO(3), where ip £ [0, 2-k), d £ [0,7r] and 
u> £ [0, 27t). We define a rotation operator on the sphere 
D(ip,d,u>) that rotates a function on the sphere, according 
to 'zyz' convention, in the sequence of u> rotation around the 
2 -axis, d rotation around the y-axis and <p rotation around 
2 -axis. A rotation of a function f(9, 0) on sphere is given by 

(V(ip,d,cj)f)(x)^f(R- 1 x), (7) 

where R is a 3 x 3 real orthogonal unitary matrix, referred 
as rotation matrix, that corresponds to the rotation operator 
D(ip,d,u>) and is given by 


where 


R = R z (ip)R y (d)R z (uj), 
cosp 

R z (p) = 


( 8 ) 


cos ip — sin ip 0 
sin ip cos ip 0 

0 0 1 


Ry(d) = 


cos d 0 sin d 
0 1 0 
— sin d 0 cos d 


( 2 ) 


is the normalization factor such that (YJ n ,Y^ = 5e tP 5 mtq , 
where S m ,q is the Kronecker delta function: 6 m ,q = 1 for 
m = q and is zero otherwise. P™(-) denotes the associated 
Legendre polynomial of degree t and order m [261. We also 
note the following relation for associated Legendre polynomial 


Here Rz{p) and R y (d) characterize individual rotations by 
p along 2 -axis and d along y-axis, respectively. 

III. Problem Lormulation 
A. Fisher-Bingham Distribution on Sphere 
Definition 1 (Lisher-Bingham Live-Parameter (LB5) Distribu¬ 
tion). The Fisher-Bingham five-parameter (FB5) distribution, 
also known as the Kent distribution, is a distribution on the 
sphere with probability density function (pdf) defined as 

0 Kp. T x+x T /3 AS 


1 


where A is a symmetric matrix of size 3x3 given by 

^ = (fjiVi vI)- 


PJfi{ cos 9) 


(3) 


(10) 
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Here fi, fji and fj -2 are the unit vectors (orthonormal set) that 
denote the mean direction (centre), major axis and minor axis 
of the distribution, respectively, k > 0 is the concentration 
parameter that quantifies the spatial concentration of FB5 
distribution around its mean and 0 < /? < k/2 is referred to 
as the ovalness parameter that is a measure of ellipticity of the 
distribution. In the term C(n , /3) denotes the normalization 
constant, which ensures ||<?|| 2 = 1 and is given by 

C(n, fi) = 2tt£ r ^ r + + 1 (^ 2r («/ 2 )~ 2r ~ 1/2 ^ + i/ 2 (*), 

( 11 ) 

where I r {') denotes the modified Bessel function of the first 
kind of order r. 

The FB5 distribution is more concentrated and more elliptic 
for larger values of k and /3, respectively. As an example, the 
FB5 distribution g(x\ K,fi, 0, A) is plotted on the sphere in 
Fig-0 for different values of parameters. 

The FB5 distribution belongs to the family of spherical 
distributions in directional statistics and is the analogue of 
the Euclidean domain bivariate normal distribution with un¬ 
constrained covariance matrix 0- ©■ 

B. 3D Spatial Fading Correlation (SFC) 

In MIMO systems, the 3D multipath channel impulse re¬ 
sponse for a signal arriving at antenna array is characterized by 
the steering vector of the antenna array. For an antenna array 
consisting of M antenna elements placed at z p £ R 3 , p = 
1,2,..., M, the steering vector, denoted by at(x), is given by 

a(x) = [a 1 (x),a 2 (x),...a L (x)], a p (x) = e lkZpX , 

( 12 ) 

where x £ R 3 denotes a unit vector pointing in the direction 
of wave propagation and k = 27t/A with A denoting the 
wavelength of the arriving signal. For any h(x) representing 
the pdf of the angles of arrival (AoA) of the multipath 
components or the unit-normalized power of a signal received 
from the direction x, the 3D SFC function between the p- 
th and the q-th antenna elements, located at z p and z q , 
respectively, with an assumption that signals arriving at the 
antenna elements are narrowband, is given by 0 

p(z p ,z q ) = / h{x) a p {x) a q (x) ds(x) 

J s 2 

= [ h(x)e ik ^- z “> a ds(x)=p(z p ~z q ), (13) 

J s 2 

which indicates that the SFC only depends on z p — z q and is, 
therefore, spatially wide-sense stationary. 

C. FB5 Distribution Based Spatial Correlation Model and 
Problem Under Consideration 

The FB5 distribution offers the flexibility, due to its direc¬ 
tional nature, to model the distribution of normalized power or 
AoA of the multipath components for most practical scenarios. 
Utilizing this capability of FB5 distribution, a 3D spatial 
correlation model for the mixture of FB5 distributions defining 
the AoA distribution has been developed in |9j. Here the 


mixture refers to the positive linear sum of a number of FB5 
distributions, each with different parameters. We defer the 
formulation of FB5 based correlation model until Section m 
Although the correlation model proposed in |9) is general in a 
sense that the SFC function can be computed for any arbitrary 
antenna geometry, the formulation of SFC function involves 
the computation of integrals that can only be carried out using 
numerical integration techniques as we highlighted earlier. 
If the spherical harmonic expansion of the FB5 distribution 
is given in closed-form, the SFC function can be computed 
analytically following the approach used in 0- 

In this paper, we derive an exact expression to compute 
the spherical harmonic expansion of the FB5 distribution. 
Using the spherical harmonic expansion of FB5 distribution, 
we also formulate exact SFC function for the mixture of FB5 
distributions, defining the distribution of AoA of the multipath 
components. We also analyse the computational considerations 
involved in the evaluation of spherical harmonic expansion of 
FB5 and distribution and SFC function. 


IV. Spherical Harmonic Expansion of the FB5 
Distribution 


In this section, we derive an analytic expression for the 
computation of spherical harmonic coefficients of the FB5 
distribution. For convenience, we determine the spherical har¬ 
monic coefficients of the FB5 distribution with mean (centre) 
//° located on 2 -axis and major axis rf( and minor axis ffj 
aligned along x-axis and y-axis respectively, that is. 


= [0 0 if rfi = [1 0 0] T rg = [01 0] T (14) 

With the centre, major and minor axes as given in (li¬ 
the FB5 distribution, referred to as the standard Fisher- 
Bingham (FB) distribution, has the pdf given by 


f(x;K,0) = 


C{n,p) 


0 k , cos 6-\-(3 sin 2 0 cos 20 


(15) 


which is related to the FB5 distribution pdf g(x; k, A) 
given in d through the rotation operator T>(tp, d, u>) as 


g{x;n,/3,fi, A) = (V(ip,i3,u)f)(x;K,l3) = f(R 1 x;k,/3), 

(! 6 ) 

where the rotation matrix is related to fi, i)\ . 772 (parameters 
of FB5 distribution) as 


R= [Vufh,U], (17) 

and the Euler angles (ip, d. oj) are related to the rotation 
matrix R through (|8j. As an example, we rotate each of the 
FB5 distribution g(x; k, fi, f3, A) plotted in Fig. [I] to obtain 
the standard Fisher-Bingham distribution fix: fi), plotted 
in Fig. [2] where we have indicated the Euler angles that 
relate the FB5 distribution and the Fisher-Bingham distribution 
through ( p~6] >. 


A. Spherical Harmonic Expansion of Standard FB distribution 

Here, we derive a closed-form expression to compute the 
spherical harmonic coefficients of the standard FB distribution, 
given in 0- Later in this section, noting the relation between 
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Fig. 1: The Fisher-Bingham five-parameter (FB5) distribution g(x: k. fi, ;3. A), given in ([fijl, is plotted on the sphere for 
parameters k, /3, as indicated and: (a) fi = [0 0 1 ] T , fji = [0 1 0 ] T , minor axis rj 2 = [1 0 0] 7 . (b) fi = [1 0 0 ] T , 
771 = [0 1 0 ] T , minor axis 772 = [0 0 l] 1 . (c) fi = [0 1 0 ] T , 771 = [10 0 ] 2 , minor axis 772 = [0 0 1 ] T . 



(a) k = 25, 0 = 10, 

= (7r/2,0,0). 



(b) K = 100 , 0 = 10 , 

(<£,$, w) = (7 t/2, 7t/2, 0). 



(c) k = 100, 0 = 49, 

= (7t/2, 7t/2, 7t/2). 


Fig. 2: The standard Fisher-Bingham (FB) distribution /( X] k, /)), given in ( | 1 5[ i. is plotted on the sphere for the concentration 
parameter k and ovalness parameter j3, as indicated. Each subfigure is related to the respective subfigure of Fig. [l] through the 
the relation with Euler angles indicated. 


standard FB distribution and FB5 distribution given in ( p~6] >, 
and the effect of the rotation operation on the spherical 
harmonic coefficients, we determine the coefficients of FB5 
distribution. 

The spherical harmonic coefficient, denoted by (/)"', of 
the standard Fisher-Bingham distribution given in ( fl5| can be 
expressed as 


(f)T = 


</(■;«, ^.17*)= [ 

Js 2 


f(x;K,/3)Y e m (x)ds(x). 

( 18 ) 


Since f(x;n,/3) is a real function, we only need to compute 
the spherical harmonic coefficients for positive orders 0 < 
to < l for each degree i > 0. The coefficients for the negative 
orders can be readily computed using the conjugate symmetry 
relation, noted in (| 6 |. 


To derive a closed-form expression for computing the spher¬ 


ical harmonic coefficients, we rewrite ( p~8| ) explicitly as 

Arm fir 

(f)T = rw~^~p\ / e KCOse Pr (cose) sin Odd 

p) Jo 


C(k,P) Jo 

p2ir 

X / 


e -im<t> e 0 sin- 9 cos 2 <f>^ 


n: 


/ 0 

m /*7r 


C(k,P) Jo 


e KCOS0 F m {6) P ™(cos 9) sin Odd, 


where 


F m (0)= g/3 sin 2 9 cos 2(ft—im(j> 

Jo 


which we evaluate as Q 


f ( B ) = i 2irIm / 2 ^ sin2 ^ 

m \0 to, G 1, 3, 5, .... 


1 Using Mathematica. 
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By expanding the modified Bessel function / m / 2 (/3sin 2 9) as 

/ o\2t+m/2 / , n\At+m 

/„ /2 (^*)=x:(f) (. 9 , 


£—0 

the exponential e Kcos 0 as (8 

K, COS 6 


I - oo 

W^^( 2 n + 1 )/„ + i /2 («;)P°( co s6»), (20) 

' K n=0 


and using the relation between associated Legendre poly¬ 
nomial and Wigner-d function given in ([3]), along with the 
following expansion of Wigner-d function in terms of complex 
exponentials 


l 

dL,nW = i n ~ m E <r>/2) < n (^r/2) (21) 

u=—£ 


Once the Euler angles {ip, 0, oj) are extracted from the 
parameters using (]24 [i-(|26[i, the spherical harmonic 

coefficients of the FB5 distribution k,A) given in |9]) 

can be computed following the effect of the rotation operation 
on the spherical harmonic coefficients as (26) 

(9)7 = <«/(•; M,«, A), YD = 0, w)f, YD 

t 

= 'E De rn,m'M“)(f)7\ (27) 

m' =—£ 

where (f)" 1 is the spherical harmonic coefficient of degree 
£ and order ml of the standard Fisher-Bingham distribution 
f(x). In ( [27] ), D e m m ,(ip,0,uj) denotes the Wigner-D function 
of degree £ and orders |m|, \m'\ < l and is given by |26 | 

D e m>rn ,(v,0, u ) = (28) 


we obtain a closed-form expression for the spherical harmonic 
coefficient in (fl 8 | as 


(f)7 = 


C(re, /3) V 2re 


2£ + l 


E!(2 n + Y)I n+ i/ 2 (n) x 


n —0 
oo 


E W,oW2)) E 


( / 3 / 2 ) 2 t + m /2 


£—0 


r(* + i)r(* +m/2 + 1) 


E d «',o( 7 r / 2 ) C C,m( 7r / 2 ) G ( 4i + TO + 1 , U + «'), 

u' ——l 

( 22 ) 

where T(-) denotes the Gamma function and we have used the 
following identity |27] Sec. 3.892] 


G(p,q)= [ (sin 6) p e lqe d9, 

Jo 


Tre ^/ 2 r(p + 2) 


2p( P + i)r(e±|± 2 )r(£=|± 2 )' 


(23) 


B. Spherical Harmonic Expansion of FB5 distribution 

We use the relation between standard FB distribution and 
FB5 distribution given in © to determine the spherical 
harmonic expansion of FB5 distribution. The Euler angles 
(p,0,ut) in ©, which characterize the relation between the 
standard FB distribution and the FB5 distribution, can be 
obtained from the rotation matrix R formulated in © in 
terms of the parameters of FB5 distribution. By comparing 
© and ( | 17[ ), 0 is given by 

0 = cos _ 1 (i? 3j3 ), (24) 


where R a 5 denotes the entry at a-th row and 6 -th column 
of the matrix R given in ( fl7j ). Similarly ip £ [0,27r) and 
u> £ [0, 27t) can be found by a four-quadrant search satisfying 


sin ip = 


sinw = 


7?2,; 


y/1 - (R 3 , 3 ) 2 
7 ( 3,2 


, cos ip = 


R 


1,3 


cos a; = 


yi - (^ 3 , 3 ) 2 

—R 3.1 

C 1 - (^3, 3 ) 2 


, (25) 

, (26) 


respectively. 


C. Computational Considerations 

Here we discuss the computation of Wigner-d functions, at 
a fixed argument of n/2, which are essentially required for 
the computation of spherical harmonic expansion of standard 
FB or FB5 distribution using the proposed formulation ( |22] i. 
Another computational consideration that is addressed here is 
the evaluation of infinite summations over t and n involved in 
the computation of ( |22| . 

1) Computation of Wigner-d functions: For the computation 
of spherical harmonic coefficients of the standard FB distri¬ 
bution using ^22) , we are required to compute the Wigner- 
d functions at a fixed argument of tv/2, that is d^ rra (7t/2), 
for each |it|, \m\ < £. Let Df denote the matrix of size 
(2£ + 1) x {21 + 1) with entries d l u m (itl2£) for |u|, |m| < L 
The matrix can be computed for each £=1,2,..., using 
the relation given in |28| that recursively computes D/ from 
D f _i. 

2) Truncation over n: The infinite sum over n arises 
from ( |20| ) where an exponential function of the concentration 
measure k and the co-latitude angle 9 is expanded as an 
infinite sum of the modified Bessel functions and Legendre 
polynomials. It is well known that the Legendre polynomials 
are oscillating functions with a maximum value of one. The 
modified Bessel function I n+ i/ 2 (k) that makes up the core 
of the expansion given in ( [20] > decays quickly (and monoton- 
ically) to zero as n —> 00 for a given k. We use this feature 
of the modified Bessel function I n+ 1 / 2 (k) to truncate the 
summation over n. We plot I n+ 1/2 («) for different values of n 
and 0 < k < 100 in Fig. [3] where it is evident that the Bessel 
function quickly decays to zero. We propose to truncate the 
summation over n in ©, or equivalently in ( [22] ), at n = N 
such that /jv+i/ 2(«0 < 10“ 16 (double machine precision). We 
approximate the linear relationship between such truncation 
level and the concentration parameter k, given by 

AT=5 k + 24, (29) 

which is also indicated in Fig. [3] This truncation level is found 
to give truncation error less than the machine precision level 
for the values of concentration parameter re in the range 0 < 
re < 100 (used in practice ©’ ©)■ We note that the truncation 
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Truncation level, n 

Fig. 3: Surface plot of modified Bessel function I n+ i/ 2 (k) for 
different values of re and n. The region I n+ i/ 2 (k) < 10” 16 
is approximated by the region below the straight line, n = 
|re + 24. Only those values I n +x/i(n) > 10~ 16 plotted above 
the straight line account for the sum in ( f20| to make it accurate 
to the machine precision level. 



Truncation level, t 

Fig. 4: Surface plot of the term inside the summation on the 
right hand side of ( [T9| ) for m = 0 and 9 = 7r/2, that is, 
S(/3, t) = on f° r different values of the ovalness parameter 
(3 and t. The region R((3,t) < 10 -16 is approximated by the 
region below the straight line, t = ||/3 + 12. 


error becomes smaller for large values of re, indicating that 
the truncation at level less than N given in ( |29j i may also 
allow sufficiently accurate computation of spherical harmonic 
coefficients. Further analysis on establishing the relationship 
between the concentration parameter re and the truncation level 
is beyond the scope of current work. 

3) Truncation over t: The expansion of modified Bessel 
functions of the first kind for a given value of /3 as given in ( fl9| ) 
introduces the infinite sum over t. Again, the decaying (non¬ 
monotonic in this case) characteristics of the term inside the 
summation on right hand side of ( p~9] > with the increase in t 
makes it possible to truncate the infinite sum given in ( fl9| that 
mainly depends on the ovalness parameter /3 with a minimal 
error. We plot the term inside the summation on the right hand 
side of ( |T9| ) for m — 0 and 9 = tt/2, that is, S(f3,t ) = 
for different values of the ovalness parameter /3 in Fig. [4] 
We again propose to truncate the summation over t at the 
truncation level T such that the terms after the truncation 
level T are less than 10“ 16 and do not have impact on the 
summation. We approximate the linear relationship between 
the ovalness parameter /3 and the truncation level T 

T = |j/3 + 12. (30) 

Since we are computing coefficients for positive orders 0 < 
m < £ for each degree £, the truncation level given by 
that is obtained for order m = 0 and 6 = 7r/2, for 
the summation in is also valid for all positive orders 
0 < m < £ and all 9 £ [0,7r] because the term inside the 
summation is maximum for m = 0 and 0 = tt/ 2 for a given 
summation variable t. 

We finally note that the truncation levels proposed here 
allow accurate computation of spherical harmonic coefficients 
to the level of double machine precision. Later in the paper, 
we validate the correctness of the derived expression for the 
computation of spherical harmonic coefficients of the standard 


FB distribution. 


V. 3D Spatial Fading Correlation for Mixture of 
Fisher-Bingham Distributions 

In this section, we derive a closed form expression to com¬ 
pute the SFC function for the spatial correlation model based 
on a mixture of FB5 distributions defining the distribution of 
AoA of multipath components Let h[x) denotes the 

pdf of the distribution defined as a positive linear sum of W 
FB5 distributions, each with different parameters, and is given 
by 


w 

^ ^ K w g (a?, k w , fi w ; fiw > A w ), (31) 

W= 1 


where K w denotes the weight of u’-th FB5 distribution with 
parameters k w , fi w , f3 w and A^,. We assume that each I\ w 
is normalized such that ||ft ||2 = 1- For the distribution with 
pdf defined in ( |3Tj ), the SFC function, given in ( fl3| has been 
formulated in |9j| and analysed for a uniform circular antenna 
array but the integrals involved in the SFC function were 
numerically computed. 

We follow the approach introduced in to derive a 
closed-form 3D SFC function using the proposed closed-form 
spherical harmonic expansion of the FB5 distribution. Using 
spherical harmonic expansion of plane waves 1291: 


OO £ 

e ikz *' £ = 47r^*^(fc||z p ||) £ Y™ (z p /\\z p \\)Yp(x), 

£—0 m ——£ 

and expanding the mixture distribution h(x) in ( |3 1 [ i. following 
Q. and employing the orthonormality of spherical harmonics, 
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we write the SFC function in ( [13] ) as m 

oo 

p(Zp-Z q ) =4TT^2i e j((k\\z p -Z q \\) X 


£=0 


Etw( I 


i=—£ 


(32) 


where 


(h)T = (h,vn 

w e 

= E E K ™ D Lm' , (33) 

w—1 m'——£ 

which is obtained by combining \21) and ( |3 1 [ i. Here, (f)T 
denotes the spherical harmonic coefficient of the standard 
FB distribution and the Euler angles (ip w ,‘d w ,oj w ) relate the 
w-th FB5 distribution g(x; K w ,p, w ,/3 W1 A W ) of the mixture 
and the standard FB distribution f(x) through ( |16| >. In the 
computation of the SFC using the proposed formulation, given 
in ( [32| , we note that the summation for I over first few terms 
yields sufficient accuracy as higher order Bessel functions 
decay rapidly to zero for points near each other in space, 
as indicated in (2TJ, We conclude this section with 

a note that the SFC function can be analytically computed 
using the proposed formulation for an arbitrary antenna array 
geometry and the distribution of AoA modelled by a mixture 
of FB5 distributions, each with different parameters. In the 
next section, we evaluate the proposed SFC function for 
uniform circular array and a 3D regular dodecahedron array 
of antenna elements. 


VI. Experimental Analysis 

We conduct numerical experiments to validate the cor¬ 
rectness of the proposed analytic expressions, formulated in 
(|22|i-(|23|) and (|32[i-(|33j). for the computation of the spherical 
harmonic coefficients of standard FB distribution and the 
SFC function for a mixture of FB5 distributions defining the 
distribution of AoA, respectively. For computing the spherical 
harmonic transform and discretization on the sphere, we em¬ 
ploy the recently developed optimal-dimensionality sampling 
scheme on the sphere gT|. Our MAT LAB based code to 
compute the spherical harmonic coefficients of the standard 
FB (or FB5) distribution and the SFC function using the results 
and/or formulations presented in this paper, is made publicly 
available. 


A. Accuracy Analysis - Spherical Harmonic Expansion of 
standard FB distribution 

In order to analyse the accuracy of the proposed analytical 
expression, for the computation of spherical harmonic coeffi¬ 
cients of the standard FB distribution, given in ( |22| ), we define 
the spatial error as 

L —1 t 

< l ) = El/(^)-EE(/)r*rM 2 , (34) 

x p £—0 m——£ 


v 10 " 5 



= 2.5, /? = 0.5 
- e- — k — 25,(3 = 10 
| - b- - « = 100,/? = 49 


0 20 40 60 80 100 120 140 


Band-limit, L 


Fig. 5: Spatial error e(L), given in ( |34[ i, between the standard 
FB distribution given in © and the reconstructed standard 
FB distribution from its coefficients up to the degree L—l. The 
convergence of the spatial error to zero (machine precision), 
as L increases, corroborates the correctness of the derived 
spherical harmonic expansion of standard FB distribution. 


which quantifies the error between the standard FB distribution 
given in © and the reconstructed standard FB distribution 
from its coefficients, computed using the proposed analytic 
expression |22|, up to the degree L—l. The summation in 
( |34] | is averaged over L 2 number of samples of the sampling 
scheme [3TJ. We plot the spatial error e(L) against band-limit 
L for different parameters of the standard FB distribution in 
Fig. a where it is evident that the spatial error converges to 
zero (machine precision) as L increases. Consequently, the 
standard FB distribution reconstructed from its coefficients 
converges to the formulation of standard FB distribution in 
spatial domain and thus validates the correctness of proposed 
analytic expression. 


B. Illustration - SFC Function 

Here, we validate the proposed closed-form expression for 
the SFC function through numerical experiments. In our anal¬ 
ysis, we consider both 2D and 3D antenna array geometries in 
the form of uniform circular array (UCA) and regular dodec¬ 
ahedron array (RDA), respectively. The antenna elements of 
M-element UCA are placed at the following spatial positions 

z p = [Rcos 2 ^,Rsin 2 ^,0] T £l 3 , (35) 

where R denotes the circular radius of the array. For the RDA, 
20 antenna array elements are positioned at the vertices of a 
regular dodecahedron inscribed in a sphere of radius R , as 
shown in Fig. [ 6 ] for R = 1. We assume that the AoA follows 
a standard Fisher-Bingham distribution. 

Using the proposed closed-form expression in ( f32] >, we 
determine the SFC between the second and third UCA an¬ 
tenna elements and plot the magnitude of the SFC function 
p(z 2 — 23 ) in Fig. [7] against the normalized radius R/A. In 
the same figure, we also plot the numerically evaluated SFC 
function, formulated in ( |T3| > and originally proposed in {9j, 
which matches with the proposed closed-form expression for 
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Fig. 6 : Regular dodecahedron array (RDA) - 20 antenna 
elements positioned at the vertices of a regular dodecahedron 
inscribed in a sphere of radius R = 1. 
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Fig. 8 : Magnitude of the SFC function p(z p — z q ) between 
antenna elements, placed at vertices shaded in black in Fig. [ 6 ] 
of 20-element RDA of radius R. 



Fig. 7: Magnitude of the SFC function p(z 2 — 23 ) for 16- 
element UCA of radius R. 


the SFC function. We emphasise that numerical evaluation 
of the integrals employs computationally intensive techniques 
to obtain sufficiently accurate results and is therefore time 
consuming. Similarly, we compute the SFC between two 
antenna elements positioned at z p and z q on RDA, which are 
indicated in Fig. [ 6 ] and plot the magnitude p(z p — z q ) Fig. [ 6 ] 
which again matches with the numerically evaluated SFC 
function given in ( fl3] l and thus corroborates the correctness 
of proposed SFC function. 

VII. Conclusions 

In this paper, the spherical harmonic expansion of the 
Fisher-Bingham distribution and a closed-form expression that 
allows the analytic computation of the spherical harmonic 
coefficients have been presented. Using the expansion of 
plane waves in spherical harmonics and the proposed closed- 
form expression for the spherical harmonic coefficients of 
the FB5 distribution, we derived an analytic formula for the 


3D SFC function between two arbitrary points in 3D-space 
when the AoA of an incident signal has the Fisher-Bingham 
distribution. Furthermore, we have validated the correctness of 
proposed closed-form expressions for the spherical harmonic 
coefficients of the Fisher-Bingham distribution and the 3D 
SFC using numerical experiments. We have focussed on the 
use of Fisher-Bingham distribution for the computation of SFC 
function. However, we believe that the proposed spherical 
harmonic expansion of the Fisher-Bingham distribution has 
a great potential of applicability in various applications for 
directional statistics and data analysis on the sphere. 
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